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The complete next-to-next-to-next-to-leading order short-distance and bound-state QCD correc¬ 
tions to T(15 l ) leptonic decay rate T(T(15') —y £ + l~) has been finished by Beneke et al. |i|. Based 
on those improvements, we present a renormalization group (RG) improved pQCD prediction for 
r(Y(lS) —y £ + £~) by applying the principle of maximum conformality (PMC). The PMC is based 
on RG-invariance and is designed to solve the pQCD renormalization scheme and scale ambigu¬ 
ities. After applying the PMC, all known-type of /3-terms at all orders, which are controlled by 
the RG-equation, are resummed to determine optimal renormalization scale for its strong running 
coupling at each order. We then achieve a more convergent pQCD series, a scheme- independent and 
more accurate pQCD prediction for T(IS') leptonic decay, i.e. Pr/isj^e+e- |pmc = 1.2701q keV, 
where the uncertainty is the squared average of the mentioned pQCD errors. This RG-improved 
pQCD prediction agrees with the experimental measurement within errors. 

PACS numbers: 13.20.-v, 13.20.He, 12.38.Bx, ll.15.Bt 


I. INTRODUCTION 

Heavy quarkonium provides an ideal platform for 
studying the non-relativistic theories, such as the non- 
relativistic Quantum Chromodynamics (NRQCD) [l| 
and the potential NRQCD (PNRQCD) @,3]. In general, 
because v 2 < v 2 and a s {rrib) < a s (m c ), the perturbative 
results for the bottomonium will be more convergent over 
the a s - and v 2 - expansion than the charmonium cases, 
where v^,c) stands for the relative velocity of constituent 
b or c quark in the bottomonium or charmonium rest 
frame. If enough bottomonium events can be generated 
at an experimental platform, we can achieve a relatively 
more definite test of those non-relativistic theories than 
the charmonium cases. 

Being an important high-energy process, the leptonic 
decay of the ground-state bottomonium T(1S) has been 
studied up to next-to-leading order (NLO) [ 3 , OJ, next- 
to-next-to-leading order (N 2 LO) 6s, 7j, and next-to-next- 
to-next-to-leading order (N 3 LO) Q- However, even by 
including the recently finished complete N 3 LO pQCD 
corrections for both the short-distance and the bound- 
state parts, the pQCD prediction for the decay rate 
Fx/is)— >e+e~ is still about 30% lower than the PDG 
value, i.e. rT(lS)^e+e-|Exp. = 1.340(18) keV 0. Even 
worse, its pQCD convergence is questionable and one 
does not know what’s the optimal behavior of the running 
coupling. It is noted that the questionable pQCD series is 
caused by using conventional scale setting, in which the 
renormalization scale is simply fixed to be ~ 3.5 GeV 
that leads to maximum decay rate and the renormaliza¬ 
tion scale uncertainty is predicted by varying it within 
the range /i r S [3,10] GeV Q. To solve such renormal¬ 
ization scale ambiguity and to improve the pQCD predic¬ 
tion, we shall use the principle of maximum conformality 


* wuxg@cqu.edu.cn 


(PMC) fl0l - fl5j to deal with T(15) leptonic decay rate up 
to N 3 LO level. 

The PMC provides a systematic procedure to set the 
optimal renormalization scale for high-energy processes 
at any order. The behavior of the running coupling is 
governed by renormalization group (RG)-equation, i.e. 
the /3-function Ml, 

oo 

/3(Us) — du s (/Tr ) /d 111 /i r — Q s ^ (/U), (1) 

2=0 

where a s = a s /47r and fi r is the renormalization scale. 
This provides the underlying principle of PMC, i.e. the 
optimal behavior of running coupling can be achieved by 
resumming all the {/3j}-terms of the process that cor¬ 
rectly determine the a s -running behavior into the cou¬ 
pling constant. Following the PMC ^-scheme, the /3- 
pattern at each perturbative order is a superposition 
of the {/3,1 -terms coming from all the lower-order a s - 
factors M- The PMC then resums the {/3j}-series ac¬ 
cording to the skeleton-like expansion that correctly re¬ 
produces the QED limit of the observable HI The re¬ 
sultant PMC scales are functions of the running coupling 
and are in general different for different orders [2lj . and 
the resultant pQCD series is thus identical to a scheme- 
independent /3 = 0 conformal series EMI. After ap¬ 
plying the PMC, the pQCD convergence can be gener¬ 
ally improved 0- One reason for such improvement lies 
in that: Being consistent with the previous treatment 
in which the /3o-series are eliminated systematically via 
the large /3o-approximation [23l - |27j . the divergent terms 


1 It is noted that there may have nj-terms (ultra-violet free and 

irrelevant to the « s -renormalization) which should be treated 
as conformal coefficients [zl ] and shall not affect our present 
PMC scale-setting. Their values may be large and may break 

the pQCD convergence in special cases. 
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(n! disappear in the PMC pQCD series due to the 

elimination of the RG-{/3i}-terms. It has been found that 
the PMC follows the RG-invariance and satisfies all the 
RG-properties In the paper, we shall show that af¬ 
ter applying the PMC, a more accurate T(1S) leptonic 
decay rate can indeed be achieved. 

The remaining parts of the paper are organized as fol¬ 
lows. In Sec.II, we will present our calculation technology 
for the Y(15) leptonic decay rate up to N 3 LO level. In 
Sec.Ill, we present numerical results. Sec.IV is reserved 
for a summary and conclusions. One appendix provide 
some computational details for PMC. 


II. CALCULATION TECHNOLOGY 

The decay rate for the channel, T(IS 1 ) -A t + t~ 
formulated as 


- x(i s)-*e+e- 


4wa 2 

9^r 




can be 


( 2 ) 


where a is the fine structure constant, m b is the 6-quark 
pole mass, and Z\ stands for the residue of the 15-wave 
two-point correlation function near (66)-threshold, which 
can be written as [29] 


Z x = |^i(0)| 2 c. 


Ex 

( , d A 

Cy 


m b 

V 3 J \ 


(3) 


where c v and d v are matching coefficients of the leading 
and sub-leading (66)-currents within the NRQCD frame¬ 
work, whose perturbative forms are 


— 1 + CfctZg, d v — 1 + dkdg, 


(4) 


fc=l 


k =l 


where a s = ct s /47t. Here |Y>i(0)| and E\ are renormalized 
wave function at the origin and binding energy of T(1S), 
which represent the bound-state contributions and also 
receive perturbative corrections from high-order heavy 
quark potentials and dynamical gluon effect, i.e. 


E 1 =E 1 


(o) 


i + E 


e k a s 


k= 1 


1^1 ( 0)| 2 = |' 01 ° ) ( O )| 2 (1 + 


(5) 


( 6 ) 


fc =i 


The LO Coulomb wavefunction at the origin and the LO 


Coulomb binding energy are given by [3CH34 
2 (m b C F a s ) 3 


( 0 ) 


87T 


1 


Ef = --m b (C F a s y 


(7) 

( 8 ) 


where C F = 4/3. As a further step, those perturbative 
coefficients e* and /,; can be separated as 


G = e? + ef 


eT, = + + 


(9) 


where G’, l nC’ and ‘us’ denote the corrections from 
the Coulomb potential, all other non-Coulomb potentials 
and ultrasoft gluon exchange, respectively. The one-loop 
and two-loop corrections for the Wilson coefficient c v 
have been given by Refs. [35l - [38j . The fermionic and the 
purely gluonic three-loop corrections to c v can be found 
in Refs. [39l - l4ll |. The one-loop correction for d v can be 
obtained from Rcf.jUj. For the bound state contribu¬ 
tions, its NLO term is from the Coulomb potential, and 
the ultrasoft correction appears first at the third order. 
Thus, we have ef = ef = 0 and ff = /f = /f = 0. 
The Coulomb, non-Coulomb and ultrasoft corrections to 
Ei and \il>i (Oil 2 have been calculated up to N 3 LO level 

in Refs.|2lli34i3- 

Up to N n LO level, one can reformulate the pQCD ap¬ 
proximate of the decay rate r-r(is)->t+£- i n a perturba¬ 
tive series as 

n 

r n = ^c i oi +3 (/i r ). (io) 

i=0 

The LO Co can be derived from Eqs. (J21HI) , and Cfi > 1) 
at each order is a combination of the coefficients Cfc, dk, 
efc and fk at different orders. There are three energy 
regions for Y(15) leptonic decay, which are characterized 
by three typical scales, i.e. the hard one Uh ~ m b , the 
soft one fx s ~ m b v b and the ultra-soft one n us ~ m b v^. 
Because v b ~ a s (m b v b ) [l|, the soft scale m b v b is usually 
replaced by m b C F a s , which is the characteristic scale of 
bottomonium and is connected to its Bohr radius via the 
relation, rsohr = 2 /(m b C F a s ). 

Practically, one can adopt any value /i/ nlt as the initial 
renormalization scale to do the renormalization, whose 
value should be large enough to ensure the pQCD calcula¬ 
tion. Under the conventional scale setting, i.e. the renor¬ 
malization scale is fixed to be fj, r = that is usually 
choose as the typical momentum of the process, the short- 
distance and bound-state corrections possess both renor¬ 
malization and factorization scale ambiguities due to the 
truncation of perturbative series. The factorization scale 
problem is another important QCD problem, es pec ially 
for the present case with several energy scales [21|. It 
has been noted that a proper choice of renormalization 
scale can lead to a smaller factorization scale depen¬ 
dence [48j]. In the paper, we shall concentrate our at¬ 
tention on solving the renormalization scale ambiguity 
and shall take the same choices for factorization scales in 
different energy regions as those suggested in the litera¬ 
ture, that is, we fix the factorization scales as: = m b , 

Us = C F a.s(us) rnh and Uus = Cpa 2 s (us)m b [29;, [H- 

SSMiiilla- 

We note that there exist logarithmic corrections such 
as the double-logarithmic In 2 a s -terms [53l, HiJ and the 
single-logarithmic In a s -terms [55|, [Hj in the perturba¬ 
tive bound-state contributions. The origin of those log¬ 
arithmic corrections is the presence of several scales in 
the threshold region. They represent a logarithm of the 
ratio of scales, e.g. a ratio of the hard scale (m b ) to 
the soft one (m b v b ) or a ratio of the soft one (m b v b ) 
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to the ultra-soft one ( rribV the resultant lnub equals 
lnct s for bound states that are approximately Coulom- 
bic, Vb oc a s j55|. These corrections are not generated by 
the renormalization group but are related to the anoma¬ 
lous dimensions of the operators in the effective Hamil¬ 
tonian Q. Thus in the following PMC treatments, the 
value of ln(a s ) is fixed and treated as conformal coeffi¬ 
cients, e.g. ln(a s ) = ln(a s (/z s )) ~ —1.1782. 

With all the known results, we are ready to do a PMC 
analysis of Y(IS') leptonic decay rate up to N 3 LO level. 
The three-loop r 3 can be written as 

r 3 = ci.o ^(/4 nit ) + (C 2 ,o + c 2 , inf )a 4 s U nit ) 

+ (c 3 ,o + C 3 ,m/ + cs^ny/a^M™*) 

+ (c 4 ,0 +C 4 .PV + C 4 , 2 ^/ +C 4 , 3 «/)a®(Mr mt )-( 11 ) 
The coefficients Cij (i > j > 0) at a certain scale can 


be read from Refs.[29l EH . EDEtl. In those references, 
the coefficients are usually given by setting the factoriza¬ 
tion scales to be equal to the renormalization scale or by 
directly setting the renormalization scale as to/,. 

PMC is a kind of /3-resummation, all RG {ft}-terms 
should be resummed to form the effective PMC scales. 
It is thus important to get the correct RG {ft}-terms 
of the process. Part of the {ft}-terms are proportional 
to the logarithmic terms as In[/xJ, mt //Zh], In[ m™VMs]> and 
In[/ij, nit /// us ], which are eliminated by specific choice of 
renormalization scale in Refs. [29l. Ell EE EE . Thus be¬ 
fore applying the PMC, one should first reconstruct all 
the coefficients with full factorization and renormaliza¬ 
tion scale dependence. This goal is achieved by using the 
scale displacement relation derived from the /3-equation 
(ED, be. the coupling 4 ) at fc t h-order can be related 
to the coupling at any other scale M 2 as 


a k M = a k M + k/3 0 In 4«* +1 (m 2 ) + k (ft In 4 + ^ 4) 

Mi V Mi 2 MI / 


0 ln ^)a^ 2 (M2) 
Mi, 


+k 


ft In 4 + 

Mi 


M 2 , 2fc + 3 a a ,_2 M 2 , (k + l)(fc + 2) a3 ^3 M 2 


-ftftift4 + 
MI 


3! 


ft 1 ft ^ 

Mi. 


a s fc+3 (M 2 ) + 0[a^ 4 (M2)]. (12) 


The derived coefficients Cij (i > j > 0) with full factor¬ 
ization and renormalization scale dependence are put in 
the Appendix. As a check of our expressions for c*j, we 
recover the Eq.(3) of Ref.[g] by taking their choices of 
Mr = Mr nlt ; M/ = Mr (m/ = M?i)Ms,Mu S ) and n f = 4, and 


by rewriting In [pi r /mb\ as In \jj, r /(mbCF<^ s {^r))] + ^ Cf + 

In a s (Mr)- 

Following the standard PMC procedures as described 
in detail in Ref. [ 03 ], we can obtain the required {ft}- 
series at each order from Eq. ED, i-e- 


ft = ri,ofl*(Mr mt ) + ( 4 * 2,0 + 3ftr 2 ,i)a 4 (M). nlt ) + (r 3) 0 + 3ftr 2 ,i + 4ftr 3i i + 6/3or 3i2 )a 3 (Mr nlt ) 

27 o o 

+(r 4 , 0 + 3ftr 2 ,i + 4ftr 3 .i + 5ftr 4i i + ^-ft/3 0 r 3! 2 + lOftrq^ + lOftrqftaf (m™*)- (13) 


The ft-coefhcients rtj (i> j > 0) can be obtained from 
the n/-coefficients Cij (i > i> 0 ) by applying basic 
PMC formulas listed in Ref.[14|. The non-conformal co¬ 
efficients Tij (j ^ 0 ) are functions of M/ mt ; while, the 
conformal coefficients r^o are independent of M/ mt - For 
convenience, we present the conformal coefficients r i 0 
with explicit factorization scale and/or initial scale de¬ 
pendence in the Appendix. 

After applying the PMC, the three-loop leptonic decay 
rate Y(IS') changes to 

r 3 = ri,oa 3 (Qi) + r 2 , 0 a 4 (Q 2 ) + r 3 , 0 a 5 s (Q 3 ) + r 4 j 0 a®(Q 4 ), 

where Qi(i = 1,2,3,4) are PMC scales at each pertur¬ 
bative order, whose expressions with explicit factoriza¬ 
tion scale and/or initial scale dependence are put in the 
Appendix. To eliminate the non-conformal ftterms, the 
renormalization scales at each perturbative order have 


been shifted from its initial value /i/ ut to the optimal 
ones Qi at different orders. The PMC scales at each or¬ 
der are determined unambiguously by resuming all the 
same type of non-conformal ftterms governed by RG- 
equation into the running coupling. The resulting pQCD 
series is identical to the one of the conformal theory with 
0 = 0 and is thus scheme independent. The PMC scales 
correctly characterize the virtuality of the propagating 
gluons and thus also allow one to determine the value of 
the effective number of flavors n /. For the present decay 
process, the number of active flavors is fixed by the num¬ 
ber of quarks in the effective theory. Since the bottom 
and the top quark have been integrated out, thus for self- 
consistency, we shall fix n/ =4 and adopt the four-flavor 
a s -running to do our discussions. Because of lacking even 
higher-order {ft}-terms, we cannot determine Q 4 , and 
we simply set Q 4 = Q 3 in the following calculation. This 
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treatment will lead to residual scale dependence, which, 
however, will be highly suppressed [ 2 lj| . 

III. NUMERICAL RESULTS 

We adopt N c = 3 for the SU (N c )-color group and 
adopt the four-loop a s -running to do the numerical anal¬ 
ysis of the Y(IS') leptonic decay rate up to three-loop 
QCD corrections. By taking a s (Mz ) = 0.1185 @, we 
obtain = 0.301 GeV. We take the fine struc¬ 

ture constant a(2m&) = 1/132.3 j57|. Using the high¬ 
est known three-loop relation between the pole mass 
and MS-running mass and taking the 6 -quark MS-mass 
rhbifhb) = 4.180 GeV we obtain the 6 -quark pole 
mass nib = 4.922 GeV@. 



FIG. 1. The decay rate F n with n = (0,1,2,3) under the 
conventional sale setting method as a function of the initial 
choice of renormalization scale /z}. nit , where T n is defined by 
Eq.Gil and stands for the decay rate with up to nth-loop 
QCD corrections. 

We first present the decay rate r-r(is)_>.f+f- with dif¬ 
ferent loop corrections in Fig. ED, in which the con¬ 
ventional scale setting method with the renormalization 
scale fi r = /4 nlt is adopted. To be self-consistent, when 
calculating T n , the (n + l)th-loop a s -running together 
with its own Aqcd value are adopted. Fig. ED agrees with 
the conventional wisdom that with the increment of loop 
corrections, the conventional scale dependence becomes 
smaller. It also indicates that the higher-order terms are 
important for an accurate pQCD prediction. 

In Fig.@, we present the initial scale dependence for 
the PMC scales Q i, Q 2 and Q 3 . Fig. ED shows that the 


2 The choice of b-quark pole mass and also (0) | 2 and E^' 1 
in final expresses ensure the correct using of PMC, since only 

those /3-terms that are pertained to the renormalization of the 
running coupling should be absorbed into the running cou¬ 
pling. Here we also do not consider the non-perturbative cor¬ 
rections/uncertainties for |?/4 0) (0)| and 



/4 nit (GeV) 


FIG. 2. The PMC scales Qi with i = (1,2,3) at each per¬ 
turbative order versus the initial renormalization scale / 4 , nlt . 
The solid, the dashed, and the dotted lines are for Qi, Q 2 , 
and Q 3 , respectively. 

PMC scales Qi are almost independent on the choice of 
initial renormalization scale /r), nlt by varying it within a 
large perturbative region such as 2 ~ 20 GeV. If setting 
/ 4 mt = mb, we find the LO PMC scale Qi — 1.31 GeV, 
the NLO PMC scale Q 2 — 2.02 GeV and the N 2 LO PMC 
scale Q 3 ~ 5.10 GeV. Those scales are different from the 
guessed value ~ 3.5 GeV that leads to maximum decay 
rate under conventional scale setting. 



LO 

NLO 

N 2 LO 

N 3 LO 

sum 

Conv. 

+0.374 

+0.125 

+0.322 

+0.061 

+0.882 

PMC 

+2.292 

-1.198 

+0.191 

-0.015 

+1.270 


TABLE I. Contributions from each order for the three-loop 
decay rate F 3 (in unit: keV) under the conventional (Conv.) 
and the PMC scale settings, respectively. p), nlt = mb- 




I<2 

I<3 

Conv. 

33.3% 

64.7% 

7.4% 

PMC 

52.3% 

17.5% 

1 . 1 % 


TABLE II. The defined K factor (K n ) for the N n LO term of 
r 3 before and after the PMC scale setting, where n = 1, 2 
and 3, respectively. /4. mt = mb- 

The non-conformal terms determine the renormaliza¬ 
tion scales at each perturbative order and the conformal 
terms as well as the resultant PMC scales accurately dis¬ 
play the magnitude of the pQCD correction at each per¬ 
turbative order. We present the contributions from each 
order for T 3 in Table ID in which the results before and 
after the PMC scale setting are presented. Under con¬ 
ventional scale setting, the N 2 LO term is about 90% of 
the LO term, and is almost three times of the NLO term, 
breaking the pQCD nature of the series. After applying 
the PMC, the pQCD convergence is improved: the mag¬ 
nitude of N 2 LO term is about 16% of the NLO term and 
the magnitude of N 3 LO term is about 8 % of the N 2 LO 
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term. This can be show more clearly by defining a K 
factor ( K n ) that equals to the magnitude of the ratio be¬ 
tween the nth-order term and the sum of all lower-order 
terms. The K factors for NLO, N 2 LO and N 3 LO terms 
are presented in Table EH 



FIG. 3. The 1% versus the initial renormalization scale /A mt 
before and after the PMC scale setting. The dashed and solid 
lines are for the conventional (Conv.) and the PMC scale 
settings, respectively. 

In Fig.®, we present the three-loop T 3 versus the 
choice of initial scale /rj, nlt , in which the results before 
and after the PMC scale setting are presented as a com¬ 
parison. Under conventional scale setting, the decay rate 
T 3 shall first increase and then decrease with the in¬ 
crement of /4 mt ; If setting /x r = ~ 3.5 GeV, we 

obtain its maximum value, which however is still lower 
than the central PDG value by about 30%. After ap¬ 
plying the PMC, the decay rate T 3 monotonously raises 
with the increment of /zj. mt , and the renormalization scale 
dependence has been greatly suppressed. By taking a 
hard enough scale such as /ij. nlt > 4 GeV, the computed 
PMC scales and the final PMC prediction for the leptonic 
T(15) decay are highly independent to its exact values. 
If taking /z). nlt = mb, we obtain 

Pt(is)^- = 1-270 keV, (14) 

which is consistent with the central PDG value within 5% 
error [g|. In Ref. Q the authors achieved a better NNLO 
prediction by including full resummation of logarithms at 
next-to-leading-logarithmic accuracy and partial contri¬ 
butions at next-to-next-to-leading logarithmic accuracy. 
The improvement of the pQCD convergence and scale 
dependence is in some sense consistent with the PMC 
prediction. This can be explained by the fact that the 
large logarithmic terms are usually accompanied by cer¬ 
tain {/3i}-terms, thus the resummation of large log-terms 
could be consistent with the PMC /3-resummation. 

For the present process, the perturbative series starts 
at a 3 -order, slight change of its argument shall result 
in large pQCD error, thus this process provides a good 
platform for testing the correct running behavior of the 


coupling constant. On the one hand, the PMC prediction 
for Y(IS') leptonic decay reads 

r T (is)^-|pMC = 1.270lo 182 ^ 0.042 ± 0-015 keXl5) 
= 1.27018:11? keV, (16) 

where the first error is the residual initial scale depen¬ 
dence for /i). nlt £ [3,10] GeV, the second error is for 
a s (Mz) Exp. = 0.1185 ± 0.0006 [j|, and the third er¬ 
ror is the estimated unknown high-order contributions. 
The errors in the second line stand for the squared aver¬ 
ages of those errors. The unknown high-order contribu¬ 
tion is predicted as ±|C3CI®|max [IB], where the symbol 
“MAX” stands for the maximum |C 3 a®| within the region 
of /r) nlt £ [3,10] GeV. This RG-improved pQCD pre¬ 
diction agrees well with the experimental measurement. 
It is noted that for the present case, even though the 
PMC scales themselves are almost flat within the region 
of aC* £ [3,10] GeV, cf. Fig.©, there is large residual 
scale dependence in com par ison to the previous PMC ex¬ 
amples, such as Refs. [58H6ll ]. Thus we need to know even 
higher-order /3-terms for this particular process so as to 
achieve accurate PMC scales and PMC predictions. 

On the other hand, the present PMC prediction on the 
Y(IS’) decay rate together with its errors can be com¬ 
pared with the prediction under the conventional scale 
setting 

r T( is)-*+/- Iconv. = 0.882t°;°80±on2 3 2 ± 0-443 ke\(17) 
= 0 . 88218:179 keV > ( 18 ) 

where the first error is initial scale dependence for /x), nlt £ 
[3,10] GeV, the second error is from a s (Mz) |e xp . un¬ 
certainty, and the third error is the estimated unknown 
higher-order contributions. The errors in the second line 
stand for the squared averages of those errors. The cen¬ 
tral decay rate is lower than the central PDG value by 
about 34%, and the much larger errors in comparison 
to the PMC prediction are caused by the large value of 
N 3 LO term at the scale 3 GeV, which are consistent with 
observation shown in Ref. fs| • 



l~Lh 

/i S 

/ins 

AT3 |conv. 

+0.004 

+0.039 

+0.002 

-0.004 

-0.033 

-0.003 

AT3|pMC 

+0.003 

-0.004 

+0.120 

-0.091 

+0.002 

-0.003 


TABLE III. The factorization uncertainties AT 3 (in units of 
keV) for the three-loop decay rate IV before and after the 
PMC scale setting, which are caused by separately varying 
fj,h, Us and n U s by ±10% of their center values, respectively. 

Let us end with a final comment on the factorization 
scale dependence. At present, we have no strict and 
systematic way to set the factorization scale, and the 
question is much more involved when there are several 
scale regions. As a reference, we present a discussion 
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- 

. Exp. 

* Conv. 

♦ PMC . 


r 3 |conv. r 3 |pMC 


FIG. 4. A comparison of T 3 together with its pQCD errors 
before and after the PMC scale setting. The theoretical errors 
are squared average of all the mentioned uncertainties. The 
PDG value, PT(is)-> e + e - |bx P . = 1.340(18) keV [{J, is included 
as a comparison. 

of factorization scale uncertainties under several simple 
choices of factorization scales, whose values before and 
after the PMC sale setting are presented in Table. IIIII 
Here, to ensure the effectiveness of the NRQCD and 
pNRQCD factorization approaches, we vary the scales 
Hh, Ms and AGs separately by ±10% of their center val¬ 
ues; e.g. when discussing the uncertainty of we take 
Hh = (1 ± 10 %)mb and fix fi s and fi us to be their cen¬ 
tral values; the uncertainties for /x s and fj, us are done via 
the same way. Table. IIIII shows that after applying the 
PMC, the factorization scale uncertainties are still there 
and the largest uncertainty is caused by the soft scale 
fi s . As shown by Eas. (IA19IA20IA2T1) . the PMC scales 
depend on the factorization scales. More explicitly, when 
setting fj, s = 90 %CFCe s (fJ-s)mb, the value of Q\ changes 
from 1.31GeV —1 1.27GeV, the value of Q 2 changes from 
2.02GeV —> 1.96GeV, and the value of Q 3 changes from 
5.10GeV -A 4.78GeV. Those are slight scale changes, 
however they shall lead to sizable contributions, since 
the decay rate starts at a^-order. 

IV. SUMMARY 

We have studied the N 3 LO short-distance and bound- 
state QCD corrections to Y(IS') leptonic decay rate of 
T(1S) -A £ + £~ by applying the PMC. A comparison of 
the three-loop r 3 together with its pQCD errors before 


and after the PMC scale setting is presented in Fig.fi]), 
where the theoretical errors are squared average of all 
the mentioned pQCD uncertainties. It shows that our 
present RG-improved pQCD prediction agrees well with 
the experimental measurement within errors. After ap¬ 
plying the PMC, the pQCD convergence of the resul¬ 
tant series is improved. Thus, the PMC does provide a 
systematic and unambiguous way to set the renormal¬ 
ization scale for any QCD processes and the accuracy 
of the pQCD prediction can be greatly improved. It is 
noted that we have not considered the non-perturbative 
corrections/uncertainties for |^i(0)| and E\, and for the 
decay rate T(T(15) -A £ + £~). Those studies shall further 
improve our present PMC predictions, which are out of 
the range of the present paper. 

Acknowledgement: This work was supported in part 
by Natural Science Foundation of China under Grant 
No.11275280, and by Fundamental Research Funds for 
the Central Universities under Grant No.CDJZR305513. 


Appendix A: The coefficients aj, the conformal 
coefficients n.o and the PMC scales for r 3 

As mentioned in the body of the text, before apply¬ 
ing the PMC scale setting, one should reconstruct all 
the coefficients with full factorization and renormaliza¬ 
tion scale dependence. In this Appendix, we first present 
the coefficients c-ij for n/-power series, and then present 
conformal coefficients r,; o and the PMC scales Qi for the 
three-loop T(15) leptonic decay rate P 3 . The full renor¬ 
malization scale and factorization scale dependence shall 
be explicitly presented. 

We take the two-loop coefficients c 3j as an example 
to explain the reconstruction procedures. As this per¬ 
turbative order, we need to deal with the two-loop QCD 
corrections to both c v and |^i(0)|. We fix the scale de¬ 
pendence for c v and |- 0 i ( 0 ) | separately, which are in dif¬ 
ferent energy regions. 

The expression of the two-loop c v (/x r = = Mfc.) 

can be found in Refs.[36|, [§7|, which involves only one 
factorization scale (ih- The expression for c v {fj, r ,^h) 
at arbitrary choice of /. i r can be derived by replacing 
o- s {Tnb) in c v (rrib,Hh) to a a (n r ) with the help of the 
scale displacement relation m- By setting /j,i = m&, 
M 2 = Mr and & = 1 in Eq. m, we obtain the required 
full factorization and renormalization scale dependence 
for c v (fi r ,[if = m?i)j which is 


Cu(Mr,M^) — 1 2 Cf~ 


5 (Mr) 


5 (Mr) 


1 LL 2 

-Cf/3o In 


mi 


(22 2 tt 2 \ (l 1 

\^-—) CFTFnh -\A CA+ Q CF 


7 t 2 C f In 
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+ 7Z CFTFn f 




13C(3) 

4 



CaC f + 



79t r 2 
36 


+ 7r 2 In 2 




,(A1) 


where C A , Cp and Tp are quadratic Casimir invari¬ 
ants [62|. For a SU(N c )-coloi group, we have Ca = N C: 
Cf = (TV 2 — 1)/2N C and Tp = 1/2. In the present case, 
rih = 1, however, we keep n/j in Eq. (EU for convenience. 
The expression of the two-loop |'0i(O )| 2 can be found 


in Ref. [43], in which the renormalization scale /x r is set to 
be the soft scale /z s . To get the expression for fi r ^ /j, s , 
we can replace a s (^ s ) by a s (fj, r ) with the help of Eq. (fl2ll : 
e.g. by setting /ri = fi s , /i 2 = /x r and k = 1 in Ea. (fl2l) . 
we obtain 


|^ 2) (0)| 2 = |^ u; (0)| 2 <J 1 + (6/3 0 T(/x r )+c^ 1 )a s ( / u r ) 




(36 + 16T(/r /l ))7r 2 C A C , p + 


178 32 L{n h ) 

9 + 3 


2/-i2 

7T G rp 


+/3g (24T 2 (/x r ) - 12L(/i fl )) + 6/3iT( Mr ) + 8/3 0 c^ 1 L( / u r 


h/>,2 


3 (Mr 


(A2) 


where /r^, and /x s stand for the hard and the soft scales, 
respectively, g and 2 are non-logarithmic parts 
of the first-order and the second-order Coulomb correc¬ 
tions [H, 501’ and for the present lS-wave bound state, 
they can be simplified as 

= 2.6229- 1.61351n/, 

c ° 2 = 1800.75 - 193.489n/ + 3.50376n^. 

The logarithmic function L{x) is defined as 

L{x) = ln[x/(mtC , pQ' s (/i s ))], (A3) 

where x = /x r , Hh or /j, s , respectively. 

As a combination, we get the required coefficients . 
Following the similar treatment, we can derive the ex¬ 
pressions for all the coefficients Cij (i > j > 0 ) with full 


scale dependence up to three-loop level from the ones at 
a particular scale presented in Refs. [29l. [jj], i43M47| . 

In using the original results of Refs.[29l. |43[ |45[, there 
are some subtleties in deriving the full scale dependence 
of the three-loop coefficient fs = / : f + c + /g s for 
I'j/’i (0)| 2 . Most of logarithmic terms for the Coulomb cor¬ 
rection /]? given in Ref. [43} are for p r = /i s , and there is 
one logarithmic term that has explicit ultrasoft scale de¬ 
pendence, which originates from the non-Abelian gluon 
“H-diagram” [ 51 ], [H} and should be treated as conformal 
coefficients. The logarithmic terms for the non-Coulomb 
correction /g C should be rewritten as ln^/m 2 and 
hi/i 2 /(mf,Cpa s (/r s )) 2 The logarithmic terms for 

the ultrasoft correction /3 s are for fi r = fi s [45|, |46j • Af¬ 
ter this clarification, we are ready to derive the full scale 
dependence of fs with the help of Eg. (THU) , which reads 


|4 3) (0)! 2 = |4 2) (0)| 2 + |4 O) (0)| 2 al^r)\n 2 C\ 
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where H i = kiCV — 1 , and for the liS-wave bound state 
the non-logarithmic part of the third-order Coulomb, 
non-Coulomb and ultrasoft corrections [H, l+l, j45| can 
be simplified as, 

c £ 3 = -39854.2 + 2005.08n/ + 19.7985n 2 + 3.61806r$, 

3 = -44754.7 - 3126.52n/, 


0,3 


= 223012. 


All the coefficients under the arbitrary choice of initial 
renormalization scale /x). nlt that are adopted in the body 
of text are in the following (in unit of GeV): 


ci , 0 = 0.0734844, 

c 2 , o = -1.37493 + 4.84997A04 nit ), 
c 2 ,i = —0.118568 — 0.293938L(/4 nit ), 

c 3 , 0 = 168.629 - 41.5323£(/4 nit ) + 213.399A 2 (/4 nit ) + 60.1699 L(/j, h ) - 106.699L(// s ) 

..in.it .. 

-34.4887 ln ^- 60.1699 ln —, 

m b m b 

..init 

C 3 = -10.7309 - 10.7761L(/tj, nit ) - 25.8665L 2 (/4. nit ) + 12.9333 L(fj, s ) + 2.09022 ln , 

m b 

c 3 ,2 = 0.257472 + 0.632361L(/tj. nit ) + 0.783834L 2 (/tj. nit ) - 0.391917L(/i g ), 

c 4 , 0 = -14311.9 - 1283.62 L(fj, h ) + 22709.0i(/rj. nit ) + 6618.69A(/t/ l )L(/xj nit ) + 832.64L 2 (/xj, nit ) 

+7824.62L 3 (/xj. nit ) + 8857.05L(/Lts) - 6618.69 L(n h )L{fi s ) - 11736.9L(/4 nit )L(/A) 

+8102.97L 2 ( A t s ) + 939.94I,(/z us ) + 518.834 ln — - 3971.21L(/xi nit ) ln — 

m b m b 

n ..init ..init 
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m b m b m b m b 
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m b 

c 4 ,i = 1458.98 - 2785.37A(^ nit ) - 401.133L(/r /l )L( A ij nit ) - 957.856L 2 (Atj. nit ) - 1422.66L 3 (Atj. nit ) 
-212.745 L(na) + 401.133L(^)L(/x s ) + 2133.99L(/t). nit )L(/ts) - 534.305L 2 (/Xs) 

ii ii ..init ..init 

+97.0847 ln — + 240.68L(/ii nit ) ln — + 400.066 ln + 275.909L(ffi nit ) ln ^— 
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..init ..init 
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—2.78696 In 2 ——, (A13) 

m b 

c 4 ,3 = 0.265871 - 1.716481/(/4 nit ) - 2.10787L 2 (/4, nit ) - 1.74185L 3 (/xj. nit ) + 0.320593L(^ S ) 

+2.61278L(^j. nit )L(At s ) - 0.261278L 2 (/x s ), (A14) 


where /j, b , fx s and /j us stand for the hard, the soft and 
the ultra-soft factorization scales, respectively. L(/4 mt ), 
L(nh)i and L(/z us ) are corresponding to taking x = 

/xj. nlt ,/xh,/z s ,/i us for L{x), respectively. 

The conformal coefficients r * q(* = 1,2,3,4) read (in 
unit of GeV), 

r lj0 = 0.0734844, (A15) 

r 2 , 0 = -3.33129, (A16) 

r 30 = 80.6951 + 60.1699A(/i ft ) - 60.1699 In -^0-17) 

m b 

r 4 ,o = 7600.61 - 1283.62L(/x h ) + 4102.54L(/x s ) 
+1634.33T 2 (/i s ) + 939.94 L{n ua ) 


+160.453L(At s ) In — + 2120.73 In — 
m b m b 

+40.1133 In 2 — - 4893.09 In — 
m b m b 

-1834.9 In 2 — + 7586.971na g (/i g ) 
m b 

+3429.11 In l n a s (^ s ) 
m b 

—2727.7In 2 a s (/i s ). (A18) 

As required, these equations show that the conformal co¬ 
efficients are free of initial scale dependence. The PMC 
scales Qi(i = 1,2,3) with full initial scale and factoriza¬ 
tion scale dependence for each perturbative order read 


i Qi 

11 (A4- nit ) 2 


11 (4 nit ) 2 


In—^1— 

(/4 nlt ) 2 
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(A19) 
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As a minor point, we have found that there are some the four-loop level, i.e. Eqs.(39b-39d) of Ref. 0 (they 

typos for the general coefficients r 4j - with j = (0,1, 2) at are correct for n = 1) should be corrected as 


r-4,2 = 
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32(n + l)(n + 2)T 3 
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